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^ . ABSTRACT 

We use a semianalytic galaxy catalogue constructed from the Millennium Simulation (MS) to 
q , study the satellites of isolated galaxies in the ACDM cosmogony. The large volume surveyed 

S-^ ' by the MS (500 3 /i _3 Mpc 3 ), together with its unprecedented numerical resolution, enable the 

C/3 ' compilation of a large sample of ~ 80,000 bright (M r < —20.5) primaries, surrounded 

by ~ 178, 000 satellites down to the faint magnitude limit (AI r = — 17) of our catalogue. 
This sample allows the characterization, with minimal statistical uncertainty, of the dynam- 
t^j- , ical properties of satellite/primary galaxy systems in a ACDM universe. The details of this 

J> ■ characterization are sensitive to the details of the modeling, such as its assumptions on galaxy 

Q> ' merging and dynamical friction timescales, but many of its general predictions should be ap- 

. plicable to hierarchical formation models such as ACDM. We find that, overall, the satellite 

' population traces the dark matter rather well: its spatial distribution and kinematics may be 

approximated by an NFW profile with a mildly anisotropic velocity distribution. Their spatial 
\Q . distribution is also mildly anisotropic, with a well-defined "anti-Holmberg" effect that reflects 

the misalignment between the major axis and angular momentum of the host halo. Our anal- 
ly , ysis also highlights a number of difficulties afflicting studies that rely on satellite velocities to 

■ constrain the primary halo mass. These arise from variations in the star formation efficiency 

and assembly history of isolated galaxies, which result in a scatter of up to ~ 2 decades in 

• i— i , halo mass at fixed primary luminosity. Our isolation criterion (primaries may only have com- 

■ panions at least 2 mag fainter within lh^ 1 Mpc) contributes somewhat to the scatter, since it 

picks not only galaxies in sparse environments, but also a number of primaries at the centre of 
"fossil" groups. We find that the abundance and luminosity function of these unusual systems 
are in reasonable agreement with the few available observational constraints. Much tighter 
halo mass-luminosity relations are found when splitting the sample by colour: red primaries 
inhabit halos more than twice as massive as those surrounding blue primaries, a difference 
that vanishes, however, when considering stellar mass instead of luminosity. The large scat- 
ter in the halo mass-luminosity relation hinders the interpretation of the velocity dispersion 
of satellites stacked according to the luminosity of the primary. We find L oc a 3 (the natu- 
ral scaling expected for ACDM) for truly-isolated primaries, i.e., systems where the central 
galaxy contributes more than 85% of the total luminosity within its virial radius. Less strict 
primary selection, however, leads to substantial modification of the scaling relation: blindly 
stacking satellites of all primaries results in a much shallower L-a relation that is only poorly 
approximated by a power law. 

Key words: galaxies: haloes - galaxies: kinematics and dynamics - dark matter - methods: 
statistical 



1 INTRODUCTION 
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Satellite galaxies may be thought of as the visible fossil relics of 
hierarchical galaxy formation, where the mass of a galaxy is envi- 
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sioned to be assembled in a sequence of accretion events. As surviv- 
ing witnesses of the accretion process, satellites bear an invaluable 
record of the assembly history of the primary galaxy they orbit, 
and provide at the same time prime information about the mass and 
extent of the dark matter halo they inhabit. 

Satellites are a particularly valuable tool for studying the 
outer regions of dark matter halos, where few other tracers ex- 
ist that can provide effective constraints. This is especially true 
in the case of the Local Group, where the dynamics of the outer 
satellites has played an important role in mass estimates of the 
Milky Way and M31 (Eittle & Tremainell987l : IZaritskv et alj 19891; 
Kochanekl fl996l: I Wilkinson & Evans! 1 19991 : lEvans & Wilkinsor] 
200G; Battaglia et al. 2005). Local Group satellites can be found 



even if they are extremely faint; Draco, for example, is ~ 8 x 10 
times fainter than its primary, the Milky Way, and some of the 



ties rivaling that of ordinary star clusters ( 


Zucker et al]|2004 


,[2006] 


Willman et alj 120051; iMartin et al.1 12006 


; Belokurov et al. 


poo3 


20071; llrwin et alj |2007|; iMaiewski et al. 


2007; llbata et al 


2007 


Chapman etai]|2007). 



On the other hand, extragalactic satellite studies have been tra- 
ditionally limited by the scarcity of satellite/primary systems easily 
accessible to observation. Partly as a result of the strict isolation 
criteria that are imposed on the primaries in order to minimize in- 
terlopers and to avoid complications that may arise from having 
more than one dominant object, it is rare for primaries to have more 
than a few satellites bright enough for spectroscopic confirmation 
and follow up. 

The scarcity of satellites surrounding one given primary has 
prompted the adoption of stacking techniques in order to overcome 
small number statistics. For example, the satellites of all primaries 
of given luminosity, L, might be combined to yield estimates of 
the average, rather than individual, halo properties as a function 
of L. These techniques are clearly vulnerable to the presence of 
luminosity-dependent biases in the satellite distribution and of sys- 
tematic trends between halo mass and primary luminosity that may 
be difficult to detect and evaluate in observational datasets. 

Such techniques, nevertheless, seem appropriate to address a 
number of important issues in galaxy formation studies. How does 
the average dark halo mass depend on the luminosity of the pri- 
mary? What is the mass profile of the dark halo and how does it 
vary with luminosity? What is the three-dimensional shape of the 
dark halo and how does it relate to the primary and to the spatial 
distribution of satellites? These are some of the questions tradition- 
ally addressed by extragalactic satellite studies, and there is a rich 
literature documenting prior attempts at exploiting the dynamics of 
satellites to constrain the mas s and extent of dark matter halos. 

The pioneering work o f IZaritskv et al] dl993l . Il997al) . for ex- 
ample, confirmed the presence of massive halos around isolated 
spirals, but also hinted at a few odd results that are difficult to rec- 
oncile with current galaxy formation models. For example, these 
authors noted that satellite velocities seemed to be independent of 
the luminosity of the primary, contrary to what might be naively 
expected from the Tully-Fisher relation. They also remarked that 
satellites seemed to be distributed anisotropicall y around the pr i- 
mary, in agreement with the early suggestion of Hol mberel d 19691) . 

These issues have been revisite d using the much larger 
datasets co mpiled by the 2 dfGRS dColless et al.1 1200 ll) and 
the SDSS jYork et al.1 l200d : IStrauss et alj |2002|) surveys, and 
broad consensus se ems to be gradually e merging. In particular, 
iMcKav et al] d2002h and lPrada et al] d2003h find that a well-defined 
trend between satellite velocities and primary luminosity does ap- 



pear when considering samples substantially larger than the ones 
considered by Zaritsky et al. Better statistics have also clarified the 
anisotropic distribution of satellites around spirals, which, contrary 
to Holmberg's suggest ion, apparently tend to avo i d the rotation axis 
of the disk terainerd [20051: lAzzaro et alj|2006l : fYang et al.ll20"06l: 
lAgustsson~& Brainerd 2007). Finally, the availability of larger sam- 
ples have also enabled studies of the satellite velocity dispersion 
profile, which may be used to probe the outer dark mass distribu- 
tion and to comp are it with cosmological N-b ody sim ulations (see 
IPrada et alj|200j also the review of lBrainerdll2004j, and van den 
Bosch et al. 2005a,b). 

The success of these studies depends crucially on identify- 
ing and understanding possibly subtle biases between dark matter 
and the satellite population. This is best accomplished through di- 
rect numerical simulations, where full 3D dynamical information 
is available and from which mock observational datasets may be 
created to assess the ability of analysis techniques to identify and 
correct for such biases. Progress, however, has been slow, mainly 
because simulations with enough dynamic range to resolve simul- 
taneously a galaxy and its s a tellites have only r e cently become 
possible (Klypinetal. 1999; Moore et al. 1999; Kravtsov ct al. 
[ 20041: iGao et alj|2004 IPiemand et al] 12004 iMaccio et alJl2O06t 
biemand et al]|2007l : iLibeskind et al]|2007l : Isales et al.ll2007bl lah. 
Further, these simulations have typically followed small volumes 
(single halos, in many cases) and many of them have focussed on 
the dark matter component only, hindering direct comparisons be- 
tween theory and observation. 

We in vestigate these issue s here using the galaxy catalogue 
created by ICroton et al] d2006h from the Millennium Simulation 
dSpringel et al]|2005l , hereafter MS for short). The catalogue is the 
result of a sophisticated semianalytic galaxy formation model ap- 
plied to an N-body simulation of unprecedented dynamic range. 
This model has been calibrated to reproduce many of the large- 
scale properties of the observed galaxy population, such as the lu- 
minosity and correlation functions, as well as their color depen- 
dence. However, no specific information about satellite systems has 
been taken into account and, therefore, the results we present here 
may be considered true theoretical "predictions" that can be con- 
trasted fruitfully with observation. 

Our principal aim is to provide a detailed characterization 
of the 3D properties of the satellite population of bright, isolated 
galaxies in the ACDM cosmogony. This characterization is in- 
tended (i) to guide the interpretation of observational datasets, (ii) 
to improve the identification of primary/satellite systems so as to 
minimize the contamination by interlopers, and (iii) to provide a 
framework within which the questions posed above may be prof- 
itably addressed. The sheer size of the catalogue, which lists ~ 10 7 
galaxies brighter than M T = —17 in the 500 3 /i~ 3 Mpc 3 volume of 
the MSB ensures that our results have negligible statistical uncer- 
tainty. 

Our plan for the paper is as follows. After a brief description 
of the MS (ii l2.lt and of the semianalytic galaxy formation model 
(§ 12. 2K we discuss our selection criteria for primaries and for satel- 
lites in § 12.31 Our main results are presented and discussed in §[3] 
We begin in § 13.11 with a discussion of the relation between halo 
mass, primary luminosity and color. Noting that our isolation crite- 
ria for primaries allows for the inclusion of "fossil groups" in our 



1 We express the present-day value of the Hubble constant as Ho = 100 h 
km/s/Mpc. Throughout this paper, masses and luminosities will assume 
h = 0.73 unless otherwise specified. 
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sample, we compare their abundance and luminosity function to ob- 
servations in § 13.21 The satellites' spatial distribution is discussed 
next (§ 13. 3K together with their kinematics (§ I3.4| |. We end the dis- 
cussion of our results by analyzing the relation between primary 
luminosity and satellite velocity dispersion (§ 13. 6t . as well as the 
presence of anisotropies (§ 13. 1\ in the satellite spatial distribution. 
We conclude with a brief summary in §|4] 



2 THE CATALOGUE 

2.1 The Millennium Simulation 

The Millennium Simulation (MS), one of the projects of the Virgo 
Consortium 0, is a cosmological N-body simulation of the ACDM 
universe that follows the evolution of more than 10 billion parti- 
cles in a box of 500 hT 1 Mpc (comoving) on a side. The parti- 
cle mass is 8.6 x 10 9 h -1 Mq, and particle-particle gravitational 
interactions are softened on scales smaller than 5 h" 1 kpc. The 
simulation adopts parameters consisten t with the WMAP1 results 
JSpergel et alj|2003l : ISeliak et alj|2005t) : fi m = 0.25, fi A = 0.75, 
h = 0.73, n — 1, and as = 0.9 (for details see Springel et 
al. 2005). Data from the simulation output at 64 times spaced 
logarithmically in expansion factor before z = 1, and at ap- 
proximately 200 Myr intervals thereafter. These data are used to 
build merger trees which encode the assembly history of each halo 
and its resolved substructure. These form the basis for the galaxy 
formation modeling. Complete data for these halo merging trees 
as well as for several galaxy formation models can be found at 
http : / / www . mpa-garching . mpg . de/millennium The Z = 
data from the Croton et al (2006) model used here can be found at 
http : //www . mpa-garching . mpg . de/galf orm/agnpaper / 



2.2 The Semianalytic Model 

The information contained in the MS is harvested using semiana- 
lytic galaxy formation models that follow, with simple, but phys- 
ically motivated, laws the formation of galaxies within the evolv- 
ing dark mat ter halos. The ac t ual pre s criptions used derive from 
the work of iKauffmann et all Jl999t) , ISpringel et all feOOlh and 
|Pe Lucia et alj J2004h . but have been reformulated in the model of 
Croton et al (2006) to take into account newer observational con- 
straints. 

In brief, the model tracks the build-up of the luminous com- 
ponent of each dark matter halo by prescribing how gas cools and 
transforms into stars, as well as how enriched gas is devolved to 
the interstellar medium in halos of various masses and at differ- 
ent stages of the hierarchy. A novel feature of the implementation 
of Croton et al (2006) is the inclusion of AGN feedback to curtail 
star formation in massive objects and to prevent the formation of 
overluminous galaxies at the center of galaxy clusters. The star for- 
mation and enrichment history of galaxies is then processed with 
standard spectrophotometric models to provide estimates of galaxy 
luminosities and colors. The main data we use here is the lumi- 
nosity in the 5 Sloan bands, ugriz, and we concentrate our analysis 
in the g and r bands. Unless otherwise specified, luminosities and 
magnitudes will refer to the r band in what follows. 

The semianalytic approach follows explicitly the evolution of 
galaxies within dark matter halos, even when they are accreted into 



(and become satellites of) larger structures. After accretion, satel- 
lite galaxies are followed consistently until their parent subhalo is 
destroyed by the tidal field of the more massive system, at which 
point the satellite orbit is subsequently identified as that of the most 
bound particle of the parent subhalo before disruption. The satellite 
is subsequently assumed to survive for some residual time consis- 
tent with dynamical friction estimates. As we shall see below, this 
careful treatment of the N-body simulation is crucial to character- 
ize the dynamics of satellites of isolated galaxies in the MS. Some 
aspects of the true evolution are nevertheless neglected; in partic- 
ular, the effects of dynamical friction on the orbit of a satellite af- 
ter its dark halo is disrupted but before it merges with the central 
galaxy of its host halo, as well as a necessarily rough estimation of 
its merging time. This neglect will affect some aspects of the re- 
sults presented below, notably the overall abundance and the radial 
distribution of satellites. 



2.3 The Primary/Satellite Galaxy Catalogue 

2.3.1 Primary galaxies 

We follow standard practice and identify a sample of primary 
galaxies from the MS galaxy catalogue through (i) an isolation cri- 
terion, imposed to ensure that a single object dominates the local 
dynamics traced by the satellites, and (ii) a brightness cutoff, im- 
posed to ensure that most primaries have a fair chance of having de- 
tectable satellites. Hereafter, we shall refer to systems brighter than 
M r = —20.5 and surrounded, within 1 hT 1 Mpc, only by com- 
panions at least 2 magnitudes fainter, as primary galaxies. (Note 
that neither the Milky Way nor M3 1 satisfy this strict isolation cri- 
terion.) 



2.3.2 Satellites 

We identify as satellites any galaxy brighter than M r = —17 
(the limiting magnitude of the catalogue) that lies within the virial 
radius of a primary. Primaries with no satellites brighter than 
M r = — 17 within their virial radius will be referred to as singles. 
Primaries with at least one satellite will be referred to as hosts. 



2.3.3 Further nomenclature 

One complication arises from the fact that, despite our isolation 
criterion, some primaries might themselves be satellites of larger 
systems. This tends to happen in the rarefied outskirts of massive 
clusters. In this case, the virial radius would refer to the cluster, 
rather than to the parent halo of the primary, and the velocities 
of nearby satellites would be contaminated by high-speed cluster 
members. To take this into account, we shall distinguish, if perti- 
nent, two classes of primaries: central and non-central. However, 
we emphasize that the fraction of non-central primaries is rather 
small, and that excluding them from the statistical analysis has no 
major influence on our main conclusions. 

A second note refers to satellites, some of which, as described 

3 We define the virial radius, T200, of a system as the radius of a sphere 
of mean density 200 times the critical density for closure, p cr i t = 



ml 



http://www.virgo.dur.ac.uk 



G ~ 277.5 h? Mq I kpc 3 . This implicitly defines the virial mass 
of a halo, M200 > as that enclosed within r'200 > ™d the virial velocity, V200 > 
as the circular velocity measured at r200 ■ Quantities characterizing a system 
will be referred to as "virial" and measured within r2uo> unless otherwise 
specified. 
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log(M !00 /M 01 ) 



Figure 1. Virial mass vs r-band luminosity relation for primary galaxies. Panel (a) shows data for all primaries, other panels are subsamples after applying a 
color cut (see labels). The total number of primaries plotted in each panel is quoted. Colors indicate the density of primaries at various locations in each panel, 
in logarithmic units; each color step corresponds to a variation of about 10 in number. The same applies to other panels, although the color coding has been 
renormalized in each case to make use of the whole color palette. Dotted, solid and dashed lines indicate loci of constant virial mass-to-light ratio, and are 
chosen to ease the comparison between panels and to outline the presence of three "sequences" in the data (see text for further discussion). 



in § 12.21 are still surrounded at the present time by their parent 
halos (hereafter WSUB satellites, for short), while others are iden- 
tified with the most bound dark matter particle of the parent sub- 
halo immediately prior to disruption (NOSUB satellites). Clearly, 
results concerning NOSUB satellites are likely to be more model- 
dependent, since the identification of the satellite as a distinct ob- 
ject depends in this case heavily on the assumptions made in the 



model about dynamical friction timescales. It is important to keep 
this distinction in mind as we try and interpret the results. 



2.3.4 Statistics 

It is interesting to assess how our primary selection criteria select 
galaxies from the general population. We find that the fraction of 
galaxies that are designated as primaries depends only weakly on 



Satellites and Fossil Groups in the Millennium Simulation 



HOST PRIMARIES 



SINGLE PRIMARIES 



0.4 
0.3 
0.3 
0.1 




0.4 



O.I 




0.4 



0.3 
0.2 
0.1 





10<lgJi<10.2 

0.11 (9.9%) 

0.11 (70. 
0.24 (19.6%):: 



1.5%):^ 



10.4<lgJi<10.6 

0.20 (1.5%)" 

0.13 (63.1%: 
0.21 (35.4%):: 



10.8<lgi<ll 

0.25 (7.2%) 

0.14 (54 
0.27 (38.6%):: 



,.2%tL 



I 



10.2<lg_L <10.4~ 

0.13 (3.7%) ~ 

0.12 (68.5%): 
0.22 (27.8%): 



10.6<lg_L<10.( 

0.21 (1.7%)~ 

0.14 (62.3%): 
0.24 (36.0%): 



lgi>ll 

L a : 

0.43 (29.1%): 

0.44 (13.4%): 
0.39 (57.5%): 



15 12 13 14 15 



0.4 
0.3 
0.2 
0.1 




0.4 



0.1 




0.4 



0.3 
0.2 
0.1 




10<lg_L<10.2 

0.17 (6.6%) 

0.09 (82 
0.23 (10.6%): 



10.4<lgi<10.6 

T 

0.29 (3.6%) 

0.11 (77 
0.16 (18.8%): 



- 



I I I I | I I I I 



10.8<lgi<ll 

0.36 (26.9%)f 

0.10 (65.1% 
0.14 (8.0%) 



10.2<lgL.<10.4" 
L 9 

0.25 (3.7%) ": 
0.10 (81.2%): 
0.19 (15.1%): 



10.6<lgL.<10.i 
T 

0.35 (6.0%K 
0.11 (75.4%) 
0.14 (18.6%) 




lgL.>ll 

0.31 (98.5%) 

(1.5%) 



lg(M 200 /M G ) 



lg(M 200 /M G 



Figure 2. Virial mass distribution in luminosity bins for host (left panel) and single (right panel) primaries, split according to color (see Figure[T}. Histograms 
have been normalized to the total number of galaxies in each luminosity bin. For hosts we have: (12, 747), (23, 585), (26, 039), (13, 817), (2, 436) and 
(454) galaxies in each of the 6 bins of increasing L, respectively. For single primaries, the corresponding numbers are : (91, 060), (58, 364), (23, 893), 
(4, 742), (361) and (66), respectively. Percentages quoted in the right column of each panel indicate the fraction of primaries in each color sequence for a 
given luminosity bin. Left columns show the rms values <r(log(M200 /Mq)) corresponding to each subsample. 



luminosity: for example, 23% of galaxies brighter than M r — —23 
are primaries, the vast majority of which (90%) are hosts; 94.6% 
of them are of the central type. Fainter galaxies (—21 < M r < 
—20.5) have a 19% chance of being classified as primaries, but the 
majority of them (85%) are actually singles and therefore will not 
contribute to satellite studies. Of these fainter primaries, 99.5% are 
classified as central. Overall, only 1% of primaries are non-central, 
and therefore our results are unlikely to be affected by the presence 
of these unusual systems. 

The final sample contains 258,321 primaries, 79,000 of 
which are hosts to at least one satellite. We find more than 178, 000 
satellites orbiting within the virial radius of the primaries, a number 
that rises to 508, 000 if all satellites within 1 ft" 1 Mpc of the pri- 
maries are considered. We list the fraction of primaries as a function 
of r-band luminosity and halo mass, along with other useful data, 
in Tables [TJ and [2] 



3 RESULTS AND DISCUSSION 

3.1 Mass-Luminosity-Color relation for primaries 

We begin by exploring the virial mass-luminosity relation for pri- 
maries in our catalogue, one of the main topics addressed by studies 
of satellite dynamics. This is shown in the top left panel of Fig- 
ure [JJ and illustrates the expected trend for brighter primaries to 
inhabit more massive halos. The most striking aspect of this figure, 
however, is the large scatter: at given luminosity, halos span over 
a decade in mass. Conversely, halos of given mass may host pri- 
maries spanning a factor of ~ 5 in luminosity. The rms deviation 



at given L is substantially smaller (see table [2} but the broad tails 
in the mass-luminosity distribution are important, and we shall see 
below. The large scatter in the mass-luminosity relation reflects the 
variety of mass assembly and star formation histories of isolated 
galaxies, and must be borne in mind when "stacking" galaxies ac- 
cording to luminosity in order to study their average halo proper- 
ties. 

Different colors in Figure [TJ code the number density of pri- 
maries at each location in the M-L plane, varying by a factor of 
~ 10 for each step in color, and decreasing from red to blue. Most 
of the galaxies are in a relatively tight "main sequence" (outlined 
by the red dots in Figure [TJ and delineated roughly by the solid 
line), but are surrounded by a broad cloud covering a large frac- 
tion of the panel. Closer inspection reveals the presence of two fur- 
ther relatively well-defined "sequences": one consisting of galaxies 
that are systematically brighter than the "main sequence", indicated 
roughly by the dotted line, and another one consisting of "underlu- 
minous" galaxies inhabiting fairly massive halos (see dashed line 
in Figure[TJ. 

We note that the "main sequence" follows roughly a constant 
virial mass-to-light ratio. This would be the natural expectation for 
halos where a similar fraction of their baryonic content has been 
transformed into stars of roughly similar mass-to-light ratio. 

The "bright sequence" (dotted line) is caused by the transient 
brightening of the luminous component resulting from a starburst 
triggered by a major merger. 

The "underluminous sequence" (dashed line), on the other 
hand, is linked to the stunted growth of the stellar component of 
central galaxies in massive halos that results from the "radio-mode" 
AGN feedback (see for details Croton et al 2006). In such systems, 
halos can increase their mass substantially while their central galax- 
ies grow only through mergers. This is indeed the way in which the 
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Figure 3. Top panel: Magnitude gap between the primary and the 10th 
brightest satellite within rzoo as a function of Lhost (expressed in units 
of h^Q Lq). The dashed line indicates the M r = —17 faint magnitude 
cutoff of our catalogue. For reference, we show also the approximate lo- 
cation of a few well-known systems, such as the Milky Way and M31 
(both of which fall outside the plot), the Virgo cluster iTrentham & Hodgkin 
|2002|) . and three " fossil" groups: RX J1340.6+40 18~l jones et al l l2000h . 
RX J1552.2+2013 iMendes de Oliveira et alj|2006l) and RX J1416.4+2315 
iCypriano et alj|2 006b. M agnitudes have been converted to the r-band us- 
ing IFukugita et al J 1 995 ) when necessary. Bottom panel: Luminosity dis- 
tribution of all red ((g — r) > 0.95) primaries (black solid histogram). The 
shaded magenta histogram indicate the subsample shown in the upper panel. 
These are the red host primaries having at least 10 satellites within T2oo ■ 



semianalytic model is able to reconcile the hierarchical growth of 
structure with the observed down-sizing of star formation in galax- 
ies and the presence of "red and dead" galaxies at early epochs 
(Croton et al 2006, De Lucia et al 2006, Bower et al 2006). 

The above interpretation suggests that the three "sequences" 
may be best appreciated by applying color cuts to the sample, since 
starburst galaxies will be bluer than average, while the opposite 
will be true for galaxies that have not formed stars recently. This 
expectation is largely borne out, as shown by the other three panels 
in FigureQ] galaxies bluer than g — r — 0.65 trace predominantly 
the "bright" sequence, whereas those redder than g — r = 0.95 
largely trace the "underluminous sequence". The "main sequence" 
is nicely traced by galaxies of intermediate, less extreme, colors. 
We hasten to add that the color discrimination is not perfect, and 
that residual evidence for the "main sequence" is clearly apparent 
in both the red and blue panels of Figure Q] 

The main conclusion from this discussion is that, at fixed lu- 
minosity, the virial mass depends strongly on color, and that a strict 
color selection criterion applied to primaries helps to tighten sub- 
stantially the scaling b etween virial mass and luminosity (see e.g., 
Klypin & Prada 2007). As we discuss below, the large scatter that 
would result from blindly stacking primaries in luminosity bins 
would most likely obscure many of the underlying trends. 

The prevalence of each of these sequences is a sensitive func- 
tion of halo mass and galaxy luminosity. This is illustrated in Fig- 
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Figure 4. Number density profile of satellites within ~ 2 virial radii, com- 
puted after scaling the position of all satellites to the virial radius of the 
primary halo, and normalized at r = r2oo- The contribution of WSUB 
and NOSUB satellites to the density profile is shown as well. Note that the 
majority of satellites near the center have had their parent halo stripped. 
The shape of the density profile is well approximated by an NFW profile 
with C200 ~ 5.6. This is slightly less concentrated than the "average" halo 
around these primaries, for which we estimate (C200) ~ 8.1. Both of these 
curves are shown in the figure, using the same normalization as for the 
satellites. Arrows indicate the radius containing half of the objects in each 
profile. 



ure|2] where each panel shows, for several luminosity bins, the dis- 
tribution of virial masses of galaxies in each of the three sequences. 
The sample is further split into host and single primaries, since, by 
definition, only the halos of host primaries are amenable to satellite 
dynamical studies. 

Blue primaries are clearly in the minority amongst host galax- 
ies, except at the brightest magnitudes, where they make a substan- 
tial (~ 29%) fraction of the brightest (L > 10 11 Lq) primaries. 
Red galaxies, on the other hand, are reasonably well represented at 
all magnitudes, making up about ~ 30-40% of all primaries fainter 
than 10 11 Lq. The "main sequence" of primaries dominates at all 
but the brightest luminosities: intriguingly, only blue and red pri- 
maries make up the population of L > 10 11 Lq hosts. These are 
clearly very unusual objects, either the result of a recent starburst 
in a relatively low mass halo, or the result of growth by mergers but 
without star formation at the center of a very massive halo. 

Single primaries show similar mass-color-luminosity trends as 
hosts, although galaxies with extremely red colors are less well rep- 
resented, and the "main sequence" is more prevalent than amongst 
hosts. Interestingly, only blue single primaries in low mass halos 
populate the brightest luminosity bin (L > 10 1 Lq): this is clearly 
the result of a major merger fueling a starburst of extreme but short- 
lived brightness. 

At given luminosity, single primaries tend to have lower virial 
masses than hosts. This is expected, since more massive halos tend 
to have more of everything, including satellites, implying a small 
bias toward more massive halos in a sample selected to have at least 
one satellite. The effect is noticeable, as witnessed by the shift in 
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Figure 5. The color dependence of the satellite number density profile. The 
thick dotted line shows the profile corresponding to all satellites, as in Fig- 
ure|4] The dashed (red) and solid (blue) curves show the contribution to the 
total profile of the reddest ({g - r) > 1.01) and bluest ((g — r) < 0.95) 
one-third of the satellites. Note that the reddest satellites are significantly 
more concentrated than the bluest ones. Arrows indicate the radius enclos- 
ing half of the satellites in each subsample. 



the position of the green histogram in corresponding panels of the 
single and host distributions shown in Figure [2] As expected from 
the argument below, it is largest in the lowest luminosity bins (i.e., 
L < 10 10 ' 4 I/©), where it amounts to a shift by a factor of ~ 2 in 
the average mass of galaxies in the "main" sequence. 



3.2 Application to "fossil groups" 

As is clear from Figures Q] and [2] our primary selection criteria 
select a number of systems in fairly massive halos, 10 -10 Mq, 
a range usually associated with galaxy groups and poor clusters 
rather than isolated galaxies. These are, indeed, systems of galaxies 
which, by chance, have evolved a gap of at least two magnitudes 
between the brightest and second brightest galaxy, and therefore 
are included in our sample. These systems are not rare; according 
to Table [TJ the brightest galaxy of about ~ 10% of all halos in the 
range 10 13 -10 15 Mq is a primary according to our definition. 

The presence of these systems in our sample has the potential 
of biasing our satellite velocities toward large values, and there- 
fore one must be careful to take this into account in the analysis. 
We shall return to this issue below (§ 13 - 6b : but we note here that 
these are the analogs of "fossil groups", systems of galaxies which 
are unusually X-ray bright (and massive) for their optical richness. 
"Fossil" groups or clusters are defined, just like our primaries, as 
systems where a large (> 2) magnitude gap exists between the 
brightest and second-brightest galaxy. 

Although the statistics of these objects are still fairly poorly 
known, their unusual properties have attracted attention from ob- 
servers and theorists alike, and there has been speculation that 
the abundance of these s ystems might be difficul t to reproduce in 
the ACDM cosmogony ( ID'Onghia & Lakdl2004l ; iD'Onghia et al.l 



20051; iMilosavlievic et al.l 120061 ; ISommer-Larser] l200d) . Subse- 
quent work has shown, however, that the predicted abundance of 
"fossil" systems is in reasonable agreeme nt with the still rather 
uncertain observational co nstraints (see, e.g. lD'Onghia et alj2007l ; 
Ivan den Bosch et aT]| 20071) . Worries remain that the atypical lumi- 
nosity function of known 'fossil' systems may be d ifficult to recon- 
cile w ith the substructure function of ACDM halos (D'Ongh ia et al.l 
120071) . but we emphasize that only three "fossil" groups have pub- 
lished luminosity functions. This, however, should change soon, as 
available X-ray data are systematically surveyed for the presence 
of "fossil" groups, and as imaging and spectroscopic surveys pro- 
vide secure luminosity functions for a growing number of "fossil" 
groups. 

Our intention here is to characterize the abundance and lumi- 
nosity function of "fossil"-like systems in our galaxy catalogue in 
order to guide the interpretation of future observational constraints. 
As stated above, these objects are not rare, since up to one in ten 
of the most massive halos may be regarded as a "fossil" system 
according to the above definitions (see Table QJ. These systems 
are not rare amongst bright galaxies either, as shown in Table [2j 
approximately 22% of galaxies brighter than 5 x W 10 Lq live in 
a "fossil"-like environment. "fossil"-like systems are more preva- 
lent amongst fainter galaxies, which tend to inhabit lower mass ha- 
los, where large magnitude gaps between the two brightest galaxy 
members are common (e.g., the Milky Way is the "ultimate fossil 
group"). 

This gives a total number density of 8.1 x 10 -6 h 3 Mpc~ 3 
for "fossil" systems in halos exceeding 10 13 ft -1 Mq, or 8.3 x 
10 -5 h 3 Mpc -3 for fossil groups whose brightest galaxy exceeds 
5 x W 10 h~ 2 Lq. This actually exceeds the (still rather uncertain) 



observational estimates of ~ 5 x 10 



2 x 10~ b h 3 Mp c' 



Jjones et al.l 120031 ; Ivikhlinin et ail 1 19991 ; iRomer et al.l 120001) for 
groups with X-ray luminosities ^ 10 43 /i^ 2 erg s _1 . We note, how- 
ever, that only about 10 fossil groups are actually known, and that 
this may very well be an underestimate. Our preliminary conclu- 
sion is that the semianalytic approach has no difficulty accounting 
for the abundance of groups with such a strong distinction between 
brightest and second-brightest galaxy. 

A further test concerns the actual luminosity distribution of 
galaxies in the groups. In a recent paper, D'Onghia et al (2007) 
argue, on the basis of N-body simulations, that "fossil" groups 
present a challenge to the ACDM scenario, since few cold dark 
matter halos are as deficient in substructure as "fossil" groups. This 
suggestion is based on the analysis of simulations that follow solely 
the dark matter component, so it is important to check the predic- 
tions of models that actually follow the formation of the galaxies. 

We adopt a simple measure of the shape of the luminosity 
distribution of galaxies in a "fossil" group, namely the magnitude 
difference between the brightest and tenth brightest galaxy mem- 
ber, AAfio. This is shown as a function of primary luminosity 
in the top panel of Figure [3] A typical cluster with a well pop- 
ulated Sche chter-like luminosity distri bution, such as Virgo, has 
AMio ~ 2 dTrentham & Ho dgkin 2002). Individual galaxies, such 
as the Milky Way or M31, have a much larger percentage of the 
total light concentrated in a single object, and as a result their 
AMio differ markedly from Virgo; we find AMio = 9 and 11 
for the Milky Way and M31, respectively dMateolll998l) . Fossil 
groups have intermediate values of AMio; for the three groups 
with published luminosity functions, the values span the range from 
3 to ~ 5 Jjonesetall feoOO; Men des de Oliveira et al j|2006l; 



ICypriano et alj2006l) 

These values are not unusual in our catalogue of isolated pri 
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Figure 6. The average satellite luminosity (expressed in units of I/host) 
is shown as a function of radius for various primary luminosity bins. Be- 
cause of our faint magnitude cutoff (M r = —17), we consider a differ- 
ent magnitude range for satellites, depending on the brightness of the host. 
The three dot-dashed lines correspond to primaries brighter than 10 10 ' 8 Lq 
(Mr^ — 22.5). Note that when considering the full magnitude range re- 
solved in the MS, (2 < M aa t — M host < 5.5) a strong luminosity seg- 
regation is clearly present: the average brightness of satellites decreases by 
a factor of ~ 2 from the center out to the virial radius. As expected, the 
magnitude of the effect decreases as the magnitude range narrows when 
considering fainter hosts. Once this effect is taken into account, the satellite 
luminosity segregation seems to be roughly independent of host brightness. 



maries, where AMio is a strong function of the primary luminosity, 
approaching ~ 2 (the minimum possible value given our isolation 
criterion) for the brightest primaries, but increasing rapidly with 
decreasing luminosity. As shown in the top panel of Figure[3] there 
is also a large scatter in AAfio; at L ~ 10 11 Lq we find values that 
go from ~ 3 to ~ 6, spanning the range observed in fossil groups. 
It would clearly be difficult to argue on the basis of this evidence 
that there is a substantial discrepancy between ACDM predictions 
and the observations of "fossil" groups. Further data are needed to 
clarify this issue further. 



3.3 Satellite density profile 

The solid circles in Figure [4] show the number density profile of 
satellites, computed after rescaling the position of each satellite to 
the virial radius of the host and stacking the full sample. The large 
number of satellites in our catalogue makes the statistical error in 
the profile negligible; bootstrap error bars are smaller than the size 
of each symbol. 

The shape of the satellite density profile is well described 
by the NFW (Navarro, Frenk & White 1996, 1997) formula, with 
C200 = 5.6. This result is rather insensitive to halo mass or lu- 
minosity, for example, splitting the sample of primaries in two by 
halo mass results in c = 5.6 ± 0.3 and 6.2 ± 0.5 for the high and 



low mass sample, respectively. A similar split in luminosity yields 
c = 6.2 ± 0.6 and c = 5.5 ± 0.4, respectively. 

Satellites are slightly less concentrated than the host dark mat- 
ter halos. The average host dark halo concentration may be found 
by averaging the concentrations of host halos, estimated from the 
mass-concentration relation of Fausti-Neto et al (in preparation): 



log(c 20 o) = 2.1 - 0.1 log(M 2OO /(M //i)- 1 ), 



(1) 



and taking into account the lognormal dispersion <r(log(c2oo)) = 
0.10. We find an average concentration of (C200) = 8.1. 

These two NFW profiles are plotted in Figure [4] showing that 
the difference is not large. For example, the half mass radius of the 
average dark halo is 0.39 r2oo, which is very similar to the radius 
that contains half of the stacked satellites, 0.42 r2oo- We conclude 
that satellites are a relatively unbiased tracer of the dark mass dis- 
tribution within the virial radius of a halo, at least for this model of 
galaxy formation and evolution. 

In a recent paper. Ictien e t al. (2006) study the radial projected 
distribution of satellite galaxies in the SDSS. These authors find 
that the projected satellite number density profile is well fitted by a 
power law: E(i?) oc R a ,v/itha = — 1.7 ± 0.1. In order to approx- 
imately mimic the selection criteria applied by Chen et al., we have 
projected all galaxies in our catalogue within 1 h~ 1 Mpc around 
each isolated primary, considering as "satellites" those within pro- 
jected distances AR < 500 h" 1 kpc and line-of-sight velocities 
AVJos < 500 km/s (notice that from our definition of isolated 
galaxies, all "satellites" are at least two magnitude fainter than 
the primary). We find a = —1.55 ± 0.08 in the distance range 
26 < R < 500ft -1 kpc, which is consistent with the results of 
Chen et al within the quoted errors. 



3.3.1 Dependence on satellite color 

Although the satellite population as a whole traces the dark mass 
reasonably well, there is a strong dependence on satellite color. 
This is shown in Figure \5\ where the profiles of the reddest 
((9 — r ) > 1-01, red dashed line) and bluest ((g — r) < 0.95, 
solid blue line) one-third of satellites is compared with the over- 
all profile presented in Figure [4] Red satellites are clearly much 
more centrally concentrated than blue ones: half of the red sample 
is contained within 0.32 r2oo, a radius that climbs to 0.71 T2oo for 
the bluest one-third of satellites. 

This difference in concentration may be traced to the assump- 
tion in the semianalytic treatment that a satellite loses its reservoir 
of hot gas (the future fuel for star formation) once it is accreted into 
a larger structure. Thus star formation in satellites declines quickly 
after accretion: the earlier a satellite was accreted the older (and 
redder) its stellar population will be. Early-accreting satellites have 
smaller turnaround radii and are likely to orbit closer to the center 
than late accreting ones, resulting in the trend shown in Figure[5] 

A trend of similar origin is shown also by the bottom two 
curves in Figure[4] which show the contribution to the satellite den- 
sity profile from satellites that have preserved (WSUB) or lost (NO- 
SUB) their parent dark halos. Early accreting satellites are more 
likely to have been more affected by tides and to have lost their 
parent halos, leading to the spatial segregation between WSUB and 
NOSUB satellites seen here. This figure also shows clearly the im- 
portance of tracking satellites in dark matter-only simulations even 
after their parent halos have been disrupted in the tidal field of the 
primary: NOSUB satellites are crucial to the satellite profile in the 
inner regions of the primary. 
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Figure 7. The radial velocity of all satellites in our sample as a function 
of radius, both scaled to the virial values of the host halo. The solid curve 
shows the median of the distribution, while dashed lines outline the 25% 
and 75% percentiles of the distribution as a function of radius. The velocity 
dispersion decreases out to the virial radius beyond which it remains ap- 
proximately constant. A "first infall" sequence of negative radial velocities 
is clearly defined outside ~ 0.5 ?"2oo- Satellites with very high radial veloc- 
ities (sometimes exceeding ~ 3 V200) are present, especially at large radii. 
This is a result of some primaries lying in the periphery of much larger 
structures, such as galaxy clusters; the large velocities are those of cluster 
members typically unrelated to the primary. 



3.3.2 Dependence on satellite brightness 

Figure [6] explores the luminosity segregation of satellites and its 
dependence on primary luminosity. The faint magnitude cutoff of 
our catalogue implies that, in order to compare meaningfully pri- 
maries of different brightness, satellites must be selected within a 
definite magnitude range. Satellites can only have absolute magni- 
tudes between the cutoff at M r = — 17 and Mtost + 2, so that 
fainter primaries have, by construction, satellites that span a nar- 
rower magnitude range. 

The three (red) dot-dashed curves in Figure [6] indicate the ra- 
dial dependence of the average satellite luminosity (in units of the 
host's) for the brightest hosts (L r > 10 10 ' 8 L ; M r < -22.5). 
When all satellites are considered (i.e., in the magnitude range 
2 < Msat — A/i lost < 5.5; bottom curve in Figure [6} a signifi- 
cant radial trend is seen: the average satellite luminosity drops by a 
factor of ~ 2 from the center out to the virial radius. The trend is, as 
expected, more difficult to detect when a narrower range in satellite 
brightnesses is imposed (upper two curves). Such a narrower range 
is needed when considering fainter primaries, for which results are 
shown by the dashed and dotted curves in Figure [6] Once this is 
taken into account the luminosity segregation we find seems to be 
independent of the brightness of the primary. 

The luminosity segregation shown in Figure [6] is most likely 
due to the effects of dynamical friction, which operate faster on 
more massive/more luminous satellites, bringing them closer to the 



center on a shorter timescale. It should be possible to contrast this 
result against observation. However, a note of caution about this 
interpretation is in order, recalling the results of Figure [4] Since 
most satellites, especially those near the center, have lost their par- 
ent halo, their present-day location is being traced by a single dark 
matter particle, and their survival depends directly on the semian- 
alytic model assumptions about dynamical friction. As such, the 
luminosity segregation shown in Figure [6] is likely to be model- 
dependent. 



3.4 Satellite Velocities 

3.4.1 Radial velocities 

Figure UJ shows the radial velocities of all satellites in our sam- 
ple as a function of their distance to the primary, after rescaling to 
the virial quantities of the host halo. The solid line traces the me- 
dian satellite velocity in radial bins; the dashed lines the 25% and 
75% of the distribution, respectively. Note that the median velocity 
within r2oo is constant and consistent with zero, as expected from 
a relaxed population in equilibrium. 

Outside r2oo, negative velocities are more prevalent, as satel- 
lites on their first approach to the primary start to dominate. These 
"first-infall" satellites delineate the negative velocity boundary at 
all radii, forming a sequence that becomes fairly obvious out- 
side ~ 0.5r2oo in Figure [7] The velocity along this sequence de- 
creases outward, and approaches zero at r ~ 3r2oo, the approx- 
imate location of the turna round radius according to the simple 
secondary infall model (see iBertschingertll 9851 ; Iwhite et alj|l993l ; 
iNavarro & White! 19931) . 

The "first-infall" sequence is an interesting feature of the ra- 
dial velocity distribution, one whose detection may be used to yield 
a direct estimate of the mass of the host halo. This is the rationale 
of a number of studies that attempt to pin down the location of the 



turnaround radius in the outskirts of galaxy grou ps and clusters by 

:y distr ibi 

19971; lDiaferiolll999l; iGeller et alJl 19991; iBiviano & Girardill2003 



jps anc 

looking at "caustics" i n the velocity distr i bution (jD iaferio & Gellei 



Rines et alj|2003l ; iDiaferio et al]l2005h . Although some progress 
has been made on this issue, the observ ational evidence remains 
elusive and its interpretation controvers ial ( Reisenegger et al j200d ; 



Drinkwater et alj200ll ; lMahdavi et alj 2005; Mohayae e et al.l2006l ; 



Gavazzi et aU2006l) . 

FigureUJoffers a possible explanation for these findings: satel- 
lites on first approach make up a relatively small fraction of sys- 
tems populating the outskirts of the halo. We investigate this quan- 
titatively in Figure [8] where we show the satellite radial velocity 
distribution as a function of radius, x = r/r^oo, normalized to the 
virial radius. 

Satellites within the virial radius have an approximately Gaus- 
sian velocity distribution with dispersion that declines with radius. 
Interestingly, the radial velocity dispersion in the inner regions is 
comparable to the virial velocity of the halo, ay ~ 0.96 V2oo- 
This result has been used by Sales et al (2007a) to argue that the 
relatively low velocity dispersion of Galactic satellites implies a 
fairly low mass for the Milky Way halo (V200 ~ 110 km/s, see 
that paper for further details). The velocity dispersion declines out- 
ward: the best fitting Gaussian for satellites with 0.5 < x < 1 is 
ay ~ 0.7 V2oo- 

Beyond the virial radius, the radial velocity distribution be- 
comes clearly asymmetric, with an excess of satellites with negative 
radial velocities. As mentioned above, this is a result of the increas- 
ing importance of satellites on their first approach to the host halo. 
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We have chosen to quantify this with a double Gaussian fit: one of 
zero mean velocity and dispersion fit to the distribution of positive 
radial velocities, and a second one whose dispersion and mean are 
chosen so that the sum of the two Gaussians match best the whole 
distribution. The two Gaussian fits are shown with dashed lines in 
Figure[8] and we shall hereafter refer to them, for short, as the "zero 
mean velocity" component (zmv) and the "infall" component (inf). 
Values quoted in each panel of figure 8| indicate the best-fit Gaus- 
si an parameters to th e tot al distributior j. In contrast with the results 
of lPrada et al, I J2006h and lPoweil J2006TK the infall patern is insensi- 
tive to halo mass; we find a non-negligible population of first-infall 
satellites around low and high-mass halos in our catalogue. 

As expected, the mean infall velocity decreases outward, from 
V r ~ -0.85 V200 at 1 < x < 1.5 to V r ~ -O.4V200 at 
2.5 < x < 3, and consistent with a turnaround radius located 
just outside ~ 3r2oo- Interestingly, the velocity dispersion of in- 
falling satellites amounts to about 15% and 25% of the virial ve- 
locity of the host halo, about a factor of two to three colder than the 
rest of the population. The prevalence of the first-infall population 
increases outward: it makes up 15%, 25%, 36%, and 40% of all 
satellites in each of the four x > 1 bins shown in Figure(8] respec- 
tively. These results may be used to improve algorithms intended 
to detect infalling galaxies in the regions surrounding groups and 
clusters. 



3.4.2 Tangential velocities 

In order to complete our characterization of satellite velocities we 
show in Figure[9]the distribution of the tangential velocity compo- 
nents, binned by distance just as in Figure(8] The spherical compo- 
nents are measured in a reference frame where the z-axis is chosen 
to coincide with the angular momentum of the host halo, so that Vj, 
carries information about a satellite's sense of rotation relative to 
the host. 

Within the virial radius the dispersion in and V$ is sub- 
stantially lower than the radial dispersion, and there are also hints 
of significant departures from Gaussianity. There is, for example, 
an excess of satellites that co-rotate with the host (i.e., V$ > 0); 
also, the Ve distribution is platykurtic; i.e., it is more centrally 
peaked and has broader wings than a Gaussian. Such departures 
from Gaussianity are more pronounced for satellites outside the 
virial radius. 

The long- and short-dashed curves in Figure [9] show the con- 
tribution to the total velocity distribution of the "infall" and "zero- 
mean velocity" components discussed in § 13.4.11 This is done by 
assigning satellites probabilistically to each of the two components, 
according to the Gaussian decomposition shown in Figure(8] Inter- 
estingly, the Ve distribution, which is not well approximated by a 
Gaussian, is found to be the result of adding the nearly Gaussian 
distributions corresponding to the "infall" and the "zero-mean ve- 
locity" components. Note, as well, that even inside T2oo the Ve dis- 
tribution differs slightly from a Gaussian. This is most likely due to 
the presence of infalling satellites within the virial radius. Although 
they are difficult to pick up in radial velocity, there are apparently 

4 We will use V{ and tji to indicate the mean and dispersion of the i- 
component velocity distribution, where i = r,8, <j>. The skewness £ is de- 
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Figure 8. Distribution of satellite radial velocities, in bins of different dis- 
tance to the primary. We consider only central hosts here in order to mini- 
mize the contribution of velocity interlopers. The radial range considered in 
each panel is labeled by the range in x = r /V200 used. Dashed lines show 
Gaussian fits to the profiles: two Gaussians are used when the distribution 
shows a strong asymmetry between negative and positive radial velocities. 
One of the Gaussians (the zmv component) is assumed to have zero mean 
velocity and to match the distribution of positive radial velocities; the pa- 
rameters of the second one (the m/ component) are then fit so as to match the 
whole distribution. The fit parameters (mean V r , dispersion ay, skewness £ 
and kurtosis k) are given in Table|3]as well as quoted in each panel. 



enough of them to modify the Ve distribution from Gaussian to 
platykurtic. 

As may be seen from the right-hand panel of Figure [9] the V4, 
distribution of the "infall" component is approximately Gaussian 
as well, and shows only a weak excess of satellites co-rotating with 
the host. The "zero-mean velocity" component, on the other hand, 
shows a much more pronounced asymmetry and a broader disper- 
sion. It is not clear at this point what causes this difference, but 
we plan to follow it up in future work. In all cases, the dispersion in 
the "infall" component is significantly lower than in the ZMV com- 
ponent. Table [3] lists a summary of the fit parameters for all these 
velocity components. 

3.4.3 Velocity anisotropy 

The velocity dispersion of the various spherical components de- 
clines gradually with radius, as shown in Figure \\0\ The biggest 
decline is seen for the radial velocity dispersion, which drops by 
almost a factor of two from the center out to the virial radius. The 
dispersion in the other components drops with radius at a different 
rate, leading to an anisotropy profile that increases from the center 
outwards, reaches a maximum of /3 max ^0.5atr~0.2 T2oo and 
declines to become almost isotropic (3 ~ just outside the virial 
radius. 

The radial dependence of the anisotropy is a reflection of the 
increasing importance of the first-infall population at larger radii. 
Because it is on its first approach, the inf population is "stretched" 
along the radial direction and has therefore a smaller radial velocity 
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Figure 9. Distribution of the tangential spherical components of the satellite velocities, Vg and Vs, shown in various distance bins, x = r/r2oo, as labelled 
in each panel. We have split the sample in two components, an "infall" component (inf, magenta long-dashed lines) and a "zero-mean velocity" component 
(ztnv, solid black line), using the double Gaussian fits to the radial velocity distribution of Figure[8]to assign probabilistically satellites to each. Dotted curves 
are Gaussian fits to each of these distributions. Parameters (mean, dispersion cr, skewness £ and kurtosis k, are quoted in each panel and in Table|3] 



dispersion at any given radius, compared with its tangential disper- 
sion (compare, e.g., the radial and tangential dispersions for the inf 
population in Figures [8] and |9). As the prevalence of this compo- 
nent increases outwards, it brings down the radial bias characteris- 
tic of the inner regions, leading to a decline of the anisotropy in the 
outskirts of the system. In support of this interpretation, we note 
that the anisotropy of WSUB satellites differs strongly from that of 
NOSUB satellites (Figure [Tot. The former still retain their parent 
halos; are therefore more likely to have been accreted into the sys- 
tem more recently; and have velocities that are more isotropic than 
the rest. 



Figure \W\ shows that the satellite velocity dispersion at the 
virial radius is about 40% its maximum value in the inner re- 
gions. This decline is in agreement with that expected for NFW 
halos and can be successfully recovered from observational sam- 
ples when contamination from interlopers is properly accounted 
for. Although early studies suggested a nearly flat satellite veloc- 
ity dispersion profiles favouring isothermal m odels for host halos 
( Mc Kay et alj|2002t iBrainerd & Speciarj|2003h . subsequent analy- 
sis suggested that it might be due to poor rem oval of background 
and foreground interlopers. Prada et al. (2003) studied the dynamic 
of satellite galaxies in the SDSS removing interlopers from their 
samples fitting a "Gaussian + constant" function to the satellite ve- 
locity distributions. These authors find that the projected velocity 
dispersion of true satellites drops to ~ 40 — 60% its maximum 
value at a projected distance of AR ~ 300ft _1 kpc from the pri- 
mary. The analysis of satellites in the 2dFGRS catalogue also show 
a declining a profile after interloper remotion, although the mea- 
sured dr op in the veloci ty dispersion with distance is somewhat 
weaker (Brainerdl l2004ah . The satellite galaxy kinematics in the 
Millennium Simulation appears broadly consistent with these ob- 
servational results. 



3.5 Halo mass profile from satellite dynamics 

Once the spatial distribution and the kinematics of the satellite pop- 
ulation have been characterized, we may use them to constrain the 
shape of the host halo mass profile. Assuming spherical symmetry, 
equilibrium, and that satellites are massless tracers of the potential, 
Jeans' equations link the potential with the velocity dispersion and 
density profiles of the satellites. Expressed in terms of a circular 
velocity, we have: 

V?(r) = -*U^ + ^r^ + m (2) 
a In r amr 

where the terms in the right hand side may be estimated from the 
results in the preceding discussion, and are summarized in the top 
panel of Figure QT] 

The implied circular velocity profile for the average halo pop- 
ulated by the satellites is shown in the bottom panel of Figure [TT] 
This V c profile has several of the same characteristics of the NFW 
profile: it rises to a maximum and then drops near the virial radius. 
The maximum circular velocity implied is Vm ax ~ 1.25 V200, and 
occurs at r- max ~ 0.3r2oo, which corresponds to a concentration 
of C200 ~ 7. We note that this is higher than the concentration 
(~ 5.6) derived from fitting an NFW profile to the number density 
of satellites (Figure|4]l, and is closer to the average concentration of 
(c) ~ 8.1 found for the host halos in the MS (see § 13. 3t . In spite of 
these differences, it is remarkable that satellites are overall reason- 
ably good tracers of the dark mass profile. These similarities have 
also been reported in N-body/gasdynamical simulations by Sales et 
al (2007a), and augur well for studies of the mass profile of galactic 
halos based on satellite data. 
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Figure 10. Upper panel: Velocity dispersion profile of central host satellites 
in spherical coordinates, where the polar z axis is chosen to coincide with 
the direction of the host angular momentum. Lower panel: The correspond- 
ing anisotropy parameter, (3 = 1 — at 2 /2oy 2 , where at' 2 = o~g + o~\. 
We distinguish in the bottom panel the contribution of satellites that have 
preserved their parent dark matter halo (WSUB, filled triangles) and those 
that have not (NOSUB, open circles). 



3.6 Primary luminosity vs satellite velocity dispersion 

The velocity dispersion of satellites is one of the prime tools used 
to investigate the mass of dark halos surrounding isolated galaxies 
and its dependence on luminosity. In hierarchical formation scenar- 
ios like ACDM, velocities vary with mass according to the virial 
definitions discussed in the footnote to § 12.3.21 which imply that 
mass scales with velocity like M200 oc V^oo- 

Due to the small number of satellites surrounding each pri- 
mary, it is necessary to combine them in some way in order to beat 
small-number statistics. One obvious choice is to bin primaries in 
a narrow luminosity range, stack their satellites, compute their ve- 
locity dispersion, and see how the dispersion varies as a function of 
the luminosity of the primary. 

If satellite velocities are an unbiased tracer of the virial veloc- 
ity, and luminosity is a reasonable proxy for mass, then one may 
expect L oc a 3 . The former assumption appears to be borne out 
by the analysis presented above, but the latter one is afflicted by 
the large scatter in the mass-luminosity relation discussed in § 13.11 
How does this affect the scaling between primary luminosity and 
satellite velocity dispersion? 

We show this in Figure [T2j where hosts are binned according 
to their r-band luminosity, and the one-dimensional velocity disper- 
sion is computed for each bin after stacking all of their satellites. 
In each panel, the solid circles show the result of this procedure 
including all satellites of hosts satisfying the condition expressed 
in the label. The short-dashed (blue), long dashed (green) and dot- 
dashed (red) curves correspond to selecting primaries according to 
the color cuts adopted in Figure [T] The dotted lines indicate, for 
reference, the L oc <r 3 and L oc a 2 scalings. 

When considering all primaries (top-left panel), the L-a re- 
lation is poorly fit by a power law, bending from approximately 



Figure 11. Top: Terms in the right-hand side of Jeans' equation relating 
satellite dynamics and the halo mass profile (eq.|2j. The dotted line shows 
a fit to ov(r) of the form: oy/V^oo = °"0 + (x/xo) exp(— (x/xq)"), 
with x = r/r200- Best-fitting dimensionless parameters are: xo = 0.068, 
a = 3/4 and erg = 0.6. Bottom: Solid line shows the average circular ve- 
locity profile of the potential sampled by the satellites, as derived from ap- 
plying Jeans' equation. Note that the circular velocity implied by the satel- 
lite population rises to a maximum before dropping near the virial radius. 
The maximum circular velocity implied is consistent with an NFW profile 
with C200 ~ 7. This is in reasonable agreement with the average concen- 
tration ((c) Kb 8.1) of host halos inhabited by satellites, and is higher than 
the concentration derived from the satellite density profile in Figure|4] 



L oc a 2 at faint luminosities to a substantially shallower scaling for 
L > 5 x 1O 1O L0. A similar departure of the L — a relation from 
a simple p ower-law scaling has been reported and discussed pre- 
viously bv lvan den Bosch et al. IHH). This is quite different from 
the naive scaling mentioned above for hierarchical models, and is 
largely due to the large scatter in the halo mass-luminosity relation 
shown in FigureQ] The presence of some very massive halos at all 
luminosities, together with their growing importance with increas- 
ing luminosity bend the relation off the natural L oc o A scaling. At 
~ 5 x 10 10 Lq the prevalence of such massive objects increases 
and the relation becomes even shallower. 



Restricting primaries by color helps; for example, for blue 
hosts ({g — r) < 0.65) we obtain a L oc a A relation that holds over 
the whole luminosity range. However, color cuts work less well for 
redder primaries: for 0.65 < (g — r) < 0.95 hosts we recover the 
L oc <t 3 relation at faint luminosities, although for galaxies brighter 
than ~ 6 x 10 10 Lq the scaling becomes again much shallower as 
a result of the growing importance of massive halos with underlu- 
minous central galaxies. The same aplies to the primaries with the 
reddest colors, although the change in the slope occurs at fainter 
luminosities (~4x 1O 1O L0). Selecting by color alone is thus not 
enough to ensure a sample of primaries with a well-defined power- 
law scaling between luminosity and velocity dispersion. 
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Figure 12. Velocity dispersion of satellites around isolated galaxies binned 
by r-band luminosity. Dispersions are computed after stacking all satellites 
within the virial radius of the host of each galaxy in the bin. Solid connected 
circles show the results for all satellites of primaries satisfying the crite- 
rion indicated by the label in each panel. Short-dashed (blue), long-dashed 
(green) and dot-dashed (red) curves correspond to primaries selected ac- 
cording to the color cuts chosen in Figure [T] Dotted lines indicate L oc cr 3 
and L oc a 2 scalings, to ease the comparison from panel to panel. Panels 
differ in the sample of primaries used. The importance of the primary within 
its own virial radius distinguishes these samples, as given by L^ oat / L S y B , 
where L EyB is the total luminosity of galaxies within the virial radius. 
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Figure 13. Same as figure [T2l but considering stellar masses rather than r- 
band luminosities. Cuts applied on M^ ost /M sys refer also to the mass in 
stars. 



One way of eliminating undeiiuminous primaries within very 
massive halos from our sample is to consider the richness of their 
surroundings. Massive halos will typically host a large number of 
galaxies and in such systems, despite our isolation criteria, which 
ensures the dominance of the central galaxy, the other members 
may contribute a significant fraction of the total luminosity. This 
is shown in the other three panels of Figure [12] where we con- 
sider only primaries making up more than, respectively, 70%, 85%, 
and 95% of the total combined luminosity of galaxies within the 
host halo, L sys . The stricter the criteria for selecting the primaries 
the nearer the scaling is to the "natural" I oc a 3 relation. Select- 
ing truly isolated galaxies thus requires more than just imposing a 
magnitude gap, but also a color cut and a conscientious survey of 
the surroundings to weed out fossil groups and poor clusters from 
the sample which may unduly bias the L-o scaling. Alternatively, 
one may try and eliminate from the analysis the brightest galaxies, 
where the contaminat ion by "fossil" sy stems is worse. This is the 
approach adopted bv lPrada et alj d2003l) . who show that the "natu- 
ral" scaling may also be recovered in that case. 

Despite the pruning of the sample, the dependence of the satel- 
lite velocity dispersion on the color of the primary remains. For our 
strictest isolation criterion, Lhost > 0.95 L sys , the velocity disper- 
sion of satellites of red galaxies is ~ 55% higher than that of blue 
galaxies of given luminosity. This is reminiscent of the observations 
of Brainerd T. (2004a,b) in the 2dFGRS who found larger velocity 
dispersion in satellites associated to early-type (red) primaries than 
satellites of late-type (blue) hosts. 

This effect is due largely to the different stellar mass-to-light 
ratios of galaxies of different colors. Indeed, as shown in FigurefOI 
the shift in velocity dispersion between primaries of different color 
basically disappears when stellar masses are considered instead of 
r-band luminosities. Using, whenever possible, stellar mass esti- 
mates rather than luminosity in order to bin galaxies is likely to 
give more robust results. Combining this with a strict isolation cri- 
terion that evaluates not only the luminosity gap between brightest 
and second-brightest galaxy but also the richness of the surround- 
ing field appears essential in order to avoid biases and to recover 
the natural scaling expected from hierarchical structure formation 



scenarios. 



3.7 Spatial anisotropics and the Holmberg effect 

An issue that has drawn recurring attention over time is whether the 
spatial distribution of satellites around primaries has anisotropies of 
particular significance. For example, the brightest satellites around 
the Milky Way seem to align on a plane perpendicu l ar to the disk 
of the Galaxy jLvnden-Belll 1 19821; iMaiewskil 1 1994 iKroupa et all 
l2005t | Libeskind et alj 12005). and a similar result seems to ap- 
ply to at least some of the satell i tes of the Andromeda galaxy 
(iKoch & GrebeH 12004 iMetz et alj 120071) . The small number of 
satellites involved in these analyses preclude s robust conclusion s 
to be drawn on the basis of the Local Group dZentner et alj|2005l 



but it i s intriguing that both Holmberg (1969) and Zaritsky_et__ 
dl997bl) find a similar effect in their samples of satellites of spiral 
galaxies. 

The advent of large datasets, mostly from the SDSS, has al- 
lowed the issue to be revisited, and the latest work suggests that 
satellites of isolated spirals tend to distribute themselves prefer- 
entially along the direction of the disk: the opposite of the effect 
claimed by Holmberg and present in the Milky Way I Brainerd 
2005 ; lAzzaro et al]|2O06t lYang et alj|2006); I Agustsson & Brainerd 
20071) . One recent paper dSales & Lambasll2004r) argued that satel- 
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lites in the 2dfGRS actually follow Holmberg's suggestion, but this 
was in error due apparently to ambiguities in the way position an- 

les are defined in the 2dfGRS database, and has now been resolved 

Yang et alfeOOrj) . 

The anti-Holmberg effect is due to misalignment between 
the angular momentum and the triaxial structure of the host ha- 
los. A long literature has now established dark matter halos to 
be triaxial objects, with a preference for nearly prolate shapes, 
and whose angular momentum is perpendicular to the major axis 
of the halo jFrenk et alj|l988l; iBullock et al Il200ll; iBulrock j ^Oo"! 
Jing & Sutoll2002h iBailin & Steinmetzll2005l ; iHopkins et aljliooa ; 
Bett et al.l2007l) . If satellites trace the shape of the dark matter halo 
and if the spin of the central galaxy disks preserves the direction of 
the halo angular momentum, then this would explain the scarcity of 
satellites near the rotation axis of the disk. 

We examine this in Figure[l4] where the anisotropy in the spa- 
tial distribution of satellites around primaries is measured by the 
distribution of the cosine of the polar angle (measured from the an- 
gular momentum axis of the SUB FIND subhalo centred on the host 
galaxy). An isotropic distribution would be horizontal in this plot, 
and it is clear that the spatial distribution of satellites is significantly 
anisotropic. Simulated satellites clearly show an "anti-Holmberg" 
effect, aligning themselves preferentially along the plane perpen- 
dicular to the angular momentum axis of the halo. 

The effect is significant but relatively weak; satellites along 
the plane outnumber those closer to the rotation axis by roughly 
2:1. The effect depends only very weakly on the luminosity of the 
primary, as shown in Figure [14] or on the relative brightness of 
the satellites. We have also checked that the effect is essentially 
independent of the color of the satellites. Quantitatively, we show 
this in Table|4] where we report the average value of |cos(6 | ) for 
various combinations of primary/satellite luminosity/color. 

A theoretical study of these alignment issues based on an 
N-body/semi-analytic model similar to but substantially smaller 
than the Millenni um Simulation has recently been carried out by 
Kang et al. ( 2007). These authors did not require their primary sys- 
tems to be isolated, nor did they require a substantial magnitude 
difference between host galaxy and satellites. As a result it is dif- 
ficult to compare their results quantitatively with our own. Never- 
theless, there are a number of results in common between our two 
studies. They found that the observed alignments between satellite 
and central galaxies were best explained by assuming the minor 
central galaxy to align with the spin of its host subhalo (as assumed 
here) rather than assuming the major axis of the galaxy to align with 
that of its subhalo, and they found alignments of similar strength to 
those we show here. Their model gave dependences of the strength 
of the alignment signal both on the colour of the host galaxy and 
on the colour of the satellites. Such trends are weak or absent for 
the MS samples we study here. This difference is most likely due to 
the much wider range of systems included in the Kang et al study. 

The weak dependence we find of the anisotropy on galaxy 
properties such as luminosity and color seems at odds with a num- 
ber of observational studies. Some of them suggest that only satel- 
lites of red primaries are anisotropically distributed (Azzaro et al. 
2006, Yang et al. 2006, Agustsson & Brainerd 2007), and that the 
effect is enhanced when considering red companions to red pri- 
maries (Yang et al 2006). We note, however, that this result may be 
affected by the large scatter in the mass-luminosity-color relation 
discussed in § 13.11 Satellites are typically searched within a fixed 
radius (of order 500 kpc), a choice that would lead to the inclu- 
sion of a larger fraction of interlopers around blue primaries, which 
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Figure 14. Spatial anisotropy of the satellite spatial distribution, as mea- 
sured by the distribution of the cosine of the polar angle, \cos(8)\, mea- 
sured from the rotation axis of the host halo. Different colors correspond 
to different primary luminosity, as labelled. Different line types correspond 
to varying the magnitude range used to select satellites, also as labelled 
(see also Figure [6). An isotropic distribution would be horizontal in this 
plot. Satellites show a well-defined "anti-Holmberg" effect; i.e., they tend 
to populate preferentially the plane perpendicular to the angular momentum 
axis of the host halo. 



tend to inhabit halos of lower mass. This may dilute the anisotropy 
in blue subsamples, explaining the observational results. 

A similar comment applies to other potential correlations re- 
ported in the literature, such as the trend for spatial anisotropies 
to decline at large radius (Brainerd 2005; Augustsson & Brainerd 
2007) or to increase with ha lo mass (Yang et al 2006). As this pa- 
per was nearing submission, Bai lin et all fc007h argued in a recent 
preprint that the alignment of satellites may significantly depend 
on the isolation criteria applied. These authors find that satellites of 
"truly" isolated SDSS primaries do show a polar excess (the origi- 
nal "Holmberg effect"). We have checked for this in our catalogue, 
but find little dependence of our conclusions on Lhost / L sys . A de- 
tailed assessment of the correlations claimed by Bailin et al (2007) 
in light of our results is beyond the scope of this paper, but it is 
clearly an issue of interest to which we plan to return in future 
work. 



4 SUMMARY AND CONCLUSIONS 

We analyse a large sample of isolated galaxies and their satel- 
lites selected from the semianalytic galaxy catalogue constructed 
by Croton et al (2006) from the Millennium Simulation (MS). The 
large number of galaxies in the catalogue, together with the large 
volume surveyed by the MS, allow us to characterize in detail the 
3D dynamical properties of the satellite population of bright iso- 
lated galaxies. Our isolation criterion stipulates that all galaxies 
within 1 h^ 1 Mpc should be at least 2 mag fainter than the primary. 
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This criterion typically selects galaxies in sparse environments, but 
also picks systems of galaxies with peculiar luminosity gaps be- 
tween the two brightest galaxies, such as "fossil groups". 
Our main conclusions may be summarized as follows: 

• The relation between the halo mass of isolated galaxies and 
their luminosity shows very large scatter (halo masses span over 
a decade in mass at given luminosity), compromising the ability 
of studies that rely on stacking satellites of galaxies of similar lu- 
minosity to probe their dark matter halos. Selecting primaries by 
color in order to eliminate the reddest and bluest primaries helps to 
tighten the mass-luminosity relation, but still a number of "underlu- 
minous" central galaxies of massive groups remain in the sample. 
These may, however, be excluded by surveying the environment 
of the primary and rejecting those in regions of anomalously high 
richness. 

• One corollary of the above conclusion is that the relation be- 
tween primary luminosity and satellite velocity dispersion is rather 
sensitive to the primary selection criteria. Stacking all satellites of 
all primaries leads to an L-a relation that is poorly approximated by 
a power law, and much shallower that the L cc a 3 scaling expected 
for hierarchical models. Only after weeding out massive halos with 
underluminous central galaxies do we recover the expected L cc a 3 
scaling. 

• Since our isolation criterion readily selects the central galaxies 
of "fossil" groups, our analysis may be used to predict the abun- 
dance of groups with unusual gaps in the luminosity of the two 
brightest galaxies. We find that about 8 to 10% of halos exceed- 
ing 10 13 M@/h would qualify as "fossil" systems, a result that 
seems consistent with the (so far rather uncertain) observational 
constraints. We examine recent claims that the luminosity function 
of "fossil" groups may be difficult to reproduce in the ACDM cos- 
mogony but find no obvious discrepancy with observational con- 
straints for the three fossil group luminosity functions in the litera- 
ture. Further data are needed to settle this issue. 

• The density profile of satellites around primaries may be well 
approximated by an NFW profile that is slightly less concentrated 
than the average dark matter profile. This conclusion is sensitive to 
the color of the satellites; red satellites are significantly more con- 
centrated than the dark matter; the opposite is true for blue satel- 
lites. We also find evidence for luminosity segregation in the satel- 
lite population; i.e., a weak tendency for satellites near the primary 
to be brighter than those further away. 

• The velocity distribution of satellites is, like the dark mat- 
ter, dominated by radial motions within the virial radius. The 
anisotropy is maximal at intermediate radii, becoming gradually 
more isotropic near the virial radius. This is a result of the radially- 
increasing contribution of satellites on their first infall onto the pri- 
mary, a population of objects with rather small dispersion in radial 
velocity whose contribution raises the importance of tangential mo- 
tions in the outskirts of the host halo. 

• Satellites are distributed anisotropically around primaries, 
with a well-defined but relatively weak preference for the plane 
perpendicular to the angular momentum of the halo (an "anti- 
Holmberg" effect). This is consistent with the latest observational 
studies, and is a direct result of misalignment between the angular 
momentum axis and the triaxial structure of dark matter halos. 

The characterization of the satellite population of isolated, 
bright galaxies we present here has several goals: (i) to guide the 
compilation of primary/satellite systems that minimize the pres- 
ence of interlopers, (ii) to facilitate the interpretation of observa- 
tional results, and (iii) to provide predictions that may be used to 



validate the semianalytic model of galaxy formation applied to the 
Millennum Simulation, as well as, more generally, the hierarchical 
nature of galaxy assembly. Many of our results are amenable to di- 
rect confrontation with observation, and it is to be hoped that such 
a comparison will provide a number of insights into galaxy for- 
mation physics, and perhaps even some challenges to the ACDM 
paradigm. 
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Table 1. Abundance of primary galaxies in the Millennium Simulation as a function of halo mass. Columns 2 and 3 indicate the median and the rms dispersion 
of the luminosity that correspond to primaries residing in halos within each mass bin. Columns 4 and 5 show the density (n pr ; m ) and the fraction (/ pr i m ) of 
halos in the Millennium Simulation that host a primary galaxy in our catalogue respectively. The last three columns list the contribution to / pr i m in each mass 
bin after applying the color cuts shown in Figure[T] 



^200 log(L r /L e ) o-(log(L r /L )) n prim / prim / g _ r (%) / s _ r (%) / 9 _ r (%) 

[Mo/h] [median] [h 3 Mp C - 3 ] % < 0.65 < 0.95 & > 0.65 > 0.95 



10 11 


- 10 12 


10.17 


0.15 


15.04 x 10" 4 


24.7 


1.6 


21.2 


1.9 


10 12 


-10 13 


10.51 


0.18 


5.50 x 10~ 4 


72.1 


0.6 


37.0 


34.5 


10 13 


- 10 14 


10.85 


0.18 


8.37 x 10- fi 


10.1 


0.1 


0.3 


9.7 


10 14 


-10 15 


11.34 


0.16 


0.34 x 10~ 6 


8.3 


0.2 


0.6 


7.5 



Table 2. Abundance of primary galaxies in the Millennium Simulation as a function of r-band luminosity. Columns 2 and 3 indicate the median and the rms 
dispersion of the halo mass where primaries of a given luminosity reside. Columns 4 and 5 show the density (n pr ; m ) and the fraction (/ pr i m ) of galaxies in 
the Millennium Simulation that are classified as primaries in our catalogue respectively. The last three columns list the contribution to each luminosity bin 
after applying the color cuts shown in Figure [T] 
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Table 3. Moments of the satellite velocity distributions shown in Figures[8]and [9]in several distance bins. The last column show /; nt , the fraction of satellites 
in each radial shell that are in their first approach to the host. 
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Table 4. Average |cos(0)| (where 8 is the polar angle measured from the halo angular momentum axis) of all satellites for different bins in host luminosity 
and color. The last three columns show the values corresponding to three equal-number subsamples: blue satellites ((g — r) sa t < 0.95), intermediate color 
satellites (0.95 < (g - r) 8a t < 1-01) and red satellites ((g - r) sat > 1.01). 
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